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It is shown that when an few-cycle, relativistically intense, p-polarized laser pulse is obliquely 
incident on overdense plasma, the surface electrons may form ultra-thin, highly compressed layers, 
with a width of a few nanometers. These electron "nanobunches" emit synchrotron radiation co- 
herently. We calculate the one-dimensional synchrotron spectrum analytically and obtain a slowly 
decaying power-law with an exponent of 4/3 or 6/5. This is much flatter than the 8/3 power of the 
BGP (Baeva-Gordienko-Pukhov) spectrum, produced by a relativistically oscillating bulk skin layer. 
The synchrotron spectrum cut-off frequency is defined either by the electron relativistic 7-factor, 
or by the thickness of the emitting layer. In the numerically demonstrated, locally optimal case, 
the radiation is emitted in the form of a single attosecond pulse, which contains almost the entire 
energy of the full optical cycle. 



I. INTRODUCTION 

High harmonics generated by an intense laser pulse 
incident on an overdense plasma surface are a promis- 
ing bright source of short wavelength radiation with a 
number of potential applications . Because of the co- 
herent character of the process, the generated harmonics 
are locked in phase and emerge in the form of attosecond 
pulses. This could recently also be shown experimentally 
[2]. Commonly, two mechanisms that lead to harmon- 
ics generation from solid density plasma surfaces are dis- 
tinguished: coherent wake emission (CWE) [3] and the 
"relativistically oscillating mirror" (ROM) 

CWE is predominant at moderately relativistic intensi- 
ties ao ^ 1 and short, but finite plasma gradient lengths 
VN/N ~ 0.01 A. Here, oq = eA/mc^ is the relativis- 
tically normalized laser potential. In this regime the 
Brunei electrons that re-enter the plasma excite electro- 
static oscillations in the overcritical plasma regions. Be- 
cause of the strong density inhomogeneity, the electro- 
static oscillations couple back to electromagnetic modes 
and thus generate harmonics. Due to this mechanism, 
CWE spectra have a cut-off at the plasma frequency cor- 
responding to the maximum density. Further, the har- 
monics have a comparatively wide emission cone [5; . 

ROM harmonics dominate for ag ^ 1, although recent 
studies ini have shown that they can be observed even 
at lower intensities when the plasma gradient length is 
about S/N/N 0.1 A. ROM spectra have no cut-off at 
the plasma frequency, but can extend to much higher 
frequencies. Their origin is basically the relativistic non- 
linearity in the electron motion. Baeva, Gordienko and 
Pukhov (BGP) jTj have developed a theory to analyt- 
ically describe the spectrum in this highly relativistic 
regime. The BGP harmonic spectrum is a power law 
with the exponent —8/3, smoothly rolling off into an ex- 
ponential decay at some frequency scaling as This 
type of spectrum could be observed experimentally [S] . 
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So far these two mechanisms have been the widely ac- 
cepted explanations for harmonics generation at over- 
dense plasma surfaces. However, there has also been 
some evidence that the two models do not tell the whole 
story. Especially under p-polarized oblique incidence, 
with oo S> 1 some spectra could be observed in numer- 
ical simulations [9] that do not fit to the predictions of 
neither of the models. 

In this paper, we take the most general approach to 
describe the harmonics radiation of the relativistic elec- 
trons: It is described as coherent synchrotron emission 
of the one dimensional slab geometry electron distribu- 
tion. For synchrotron radiation, the optimal efficiency is 
reached, when all the radiating electrons gather in one 
extremely dense and narrow bunch. We demonstrate via 
PIC simulation, that this case can actually closely be 
achieved in relativistic laser-plasma interaction. 



II. ASSUMPTIONS AND PREDICTIONS OF 
BGP THEORY 

We start with a brief review of the BGP theory. 
The theory is one dimensional, oblique incidence can be 
treated in a Lorentz transformed frame |10| . This is ac- 
ceptable as long as the focal spot size is big compared 
to the laser wavelength [TT]. The reflection of the light 
is done by the electrons close to the surface. In general, 
the surface electrons perform highly complicated motions 
and it is a hopeless endeavour to try to describe it ana- 
lytically even in ID. Therefore we have to resort to some 
approximation. The starting point of the BGP model is 
the following, simplified boundary condition: 

E., {t - XARp(t)/c) +Er{t + XARP (t)/c) =0 (1) 

where xarp is the so called "Apparent Refiection Point" 
(ARP) position (see Eq. (18) in [7J). The whole informa- 
tion about the plasma motion has been condensed to only 
one real function XARp{t). However, we have to empha- 
size here, that Eq. (jlj is an approximation. Although it 
is not yet possible to give a fully general criterion for the 
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validity of Eq. ([T|) in terms of the initial laser and plasma 
parameters, it has been demonstrated in |12j that the 
approximation is valid if the plasma skin layer evolution 
time T is long compared to the skin length 6 in the sense 
cr S> (5 . This is the case if 5 = Ne/agNc ^ 1, or 
generally if the electron density profile is approximately 
step-like during the main interaction phase. 

Further note that Eq. ([T]) makes sense only for an ARP 
moving at a physical velocity of IccarpI < c. In this 
case it is clear from Eq. ([ij, that the reflected field is 
nothing but a phase modulation of the negative of the 
incident field. This can easily be checked inside a PIC 
simulation. Fig. [TJi shows an exemplary result. For this 
simulation, the laser was normally incident on a sharp 
edged plasma density profile. Here, Eq. ([T]) is fulfilled 
to a good approximation as can be verified from Fig. [l^: 
Apart from minor deviations, the extremal values and the 
sequence of the monotonic intervals agree, the refiected 
field is approximately a phase modulation of the incident 
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Starting from Eq. ([TJ, it is now possible to calculate 
the spectrum of the refiected radiation with only little 
knowledge of the function XARp{t). Therefore we per- 
form the Fourier transformation of Er{t) employing the 
method of stationary phase. At the stationary phase 
point, the relativistic 7-factor corresponding to the ARP 
7ARP = (l~(iniii(i;ARp/c))^)~^^^ possesses an extremely 
sharp maximum, why we also refer to this point as the 
7-spike. The high frequency spectrum does only depend 
on the surface behavior near this 7-spike. Therefore, we 
can also relax the condition for the validity of the the- 
ory: Actually, we require Eq. ([TJ only to be fulfilled in 
the neighbourhood of the 7-spike. 

The exact result of the calculation of the spectrum is 
given by Eqs. (34)-(36) in |7], in a very good approxima- 
tion it can be simplified to: 
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(2) 



wherein Ai refers to the well known Airy-function. The 
roll-off frequency LUr marks the point, where the initial 
power law decay of the spectral envelope merges into a 
more rapid exponential decay. In the region uj <^ LUr, the 
Airy function is 0(1) and one obtains the — 8/3-power 
law spectrum. For uj ^ LUr, the decay is an exponential 
one. From the integration, it also follows that ujr oc 7ar,pj 
where 7arp is taken at its maximum. The scaling with 
7arp contrast to the long known case of refiection 
at a mirror moving with constant velocity generates a 
Doppler upshift of only (1 + - /3) « 47^. 

Fig. [1]d shows the spectrum corresponding to the field 
from Fig.[l^. As expected, it agrees well with the —8/3- 
power law. 




Figure 1: Radiation in time ((a) and (c)) and spectral ((b) 
and (d)) domain for two simulations, (a) and (b) correspond 
to the BGP case: normal incidence, plasma density A^'e = 
250A'^c, sharp edged profile; (c) and (d) correspond to the 
nanobunching case: plasma density ramp oc exp(a;/(0.33 A)) 
up to a maximum density of A^e = 95 Nc (lab frame), oblique 
incidence at 63° angle (p-polarized) . Laser field amplitude 
is ao — 60 in both cases. In all frames, the reflected fleld 
is represented by a blue line. In (a) and (c), the green line 
represents the field of the incident laser and the black dashed 
lines mark the maximum field of it. In (b) and (d), the dotted 
black line represents an 8/3 power law, the red dashed line 
corresponds to the ID synchrotron spectrum Eqs. ^ and 
with ujra ~ SOOciJo and ujrf ~ 225 ujq. 



III. A NEW REGIME OF RELATIVISTIC 
HARMONICS GENERATION 

Let us now have a look at Fig. [TJ;, which shows the 
result of another PIC simulation. In this simulation the 
plasma density profile was extended over a few fractions 
of a laser wavelength, and the p-polarized laser was inci- 
dent at an angle of a = 63°. The laser pulse used was 
the same as in the first simulation. 

It is evident, that the maximum of the refiected field 
reaches out about an order of magnitude above the am- 
plitude of the incident laser. The refiected radiation can 
clearly not be obtained from the incident one just by 
phase modulation. By this we conclude, that the bound- 
ary condition Eq. ([T) fails here. Consequently, the spec- 
trum deviates from the —8/3-power law, compare Fig.[lji. 
Indeed we see that the efficiency of harmonic generation 
is much higher than estimated by the BGP calculations: 
about two orders of magnitude at the hundredth har- 
monic. Also, we can securely exclude CWE as the re- 
sponsible mechanism, since this would request a cut-off 
around uj — IOwq- From this we conclude that the radi- 
ation cannot be attributed to any of the known mecha- 
nisms. 

To get a picture of the physics behind, let us have a 
look at the motion of the plasma electrons that generate 
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Figure 3: (Color online) Formation of the nanobunch in the 
simulation corresponding to Figs.[T]:-d and[2]3. We depict the 
electron density in units of the critical density Nc, the 
transverse magnetic field component and the longitudi- 
nal electric field component in relativistically normalized 
units. 



Figure 2: (Color online) The electron density and contour 
lines (cyan, or light grey) of the emitted harmonics radiation 
for uj/iuo > 4.5, in (a) the BGP and (b) the synchrotron (or 
nanobunching) regime. The small windows inside the main 
figures show the detailed density profile at the instant of har- 
monic generation. All magnitudes are taken in the simulation 
frame. The parameters used are the same as in the other fig- 



the radiation. Figure [2] shows the evolution of the elec- 
tron density in both our sample cases. In addition to the 
density, contour lines of the spectrally filtered reflected 
radiation are plotted. These lines illustrate where the 
main part of the high frequency radiation emerges. We 
observe that in both cases the main part of the harmon- 
ics is generated at the point, when the electrons move to- 
wards the observer. This shows again that in both cases 
the radiation does not stem from CWE. For CWE har- 
monics, the radiation is generated inside the plasma, at 
the instant when the Brunei electrons re-enter the plasma 
[3j. Apart from that mutuality, the two presented cases 
appear to be very different. 



Figure [2^ corresponds to the BGP case. It can be seen 
that the density profile remains roughly step-like during 
the whole interaction process and the plasma skin layer 
radiates as a whole. This explains why the BGP theory 
works well here, as we have seen before in Fig. [TJl and b. 

However, figure ^jp looks clearly different. The den- 
sity distribution at the moment of harmonics generation 
is far from being step-like, but possesses a highly dense 
(up to ^ 10000 Nc density) and very narrow (5-like peak, 
with a width of only a few nanometers. This electron 
"nanobunch" emits synchrotron radiation coherently. 

The radiation is emitted by a highly compressed elec- 
tron bunch moving away from the plasma. However, the 
electrons first become compressed by the relativistic pon- 
deromotive force of the laser that is directed into the 
plasma, compare the blue lines in Fig. |3] During that 
phase, the longitudinal electric field component grows 
until the electrostatic force turns around the bunch, com- 
pare the green lines in Fig. [Sj Normally, the bunch will 
loose its compression in that instant, but in some cases, 
as in the one considered here, the fields and the bunch 
current match in a way that the bunch maintains or even 
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increases its compression. The final stage is depicted by 
the red Hnes in Fig. [3] 

We emphasize, that such extreme nanobunching does 
not occur in every case of p-polarized obHque incidence 
of a highly relativistic laser on an overdense plasma sur- 
face. On the contrary, it turns out that the process is 
highly sensitive to changes in the plasma density pro- 
file, laser pulse amplitude, pulse duration, angle of inci- 
dence and even the carrier envelope phase of the laser. 
For a longer pulse, we may even observe the case, that 
nanobunching is present in some optical cycles but not in 
others. The parameters in the example were selected in 
a way to demonstrate the new effect unambiguously, i.e. 
the nanobunch is well formed and emits a spectrum that 
clearly differs from the BGP one. The dependence of the 
effect on some parameters is discussed in section [Vl] 

Because of the one dimensional slab geometry, the 
spectrum is not the same as the well known synchrotron 
spectra [131 of a point particle. We now calculate the 
spectrum analytically. 



IV. SPECTRUM OF ID COHERENT 
SYNCHROTRON EMISSION (CSE) 

The radiation field generated to the left of a one- 
dimensional current distribution can in a completely gen- 
eral way be expressed as: 



Esy (t, X) 



2tt 



+ 00 



j {t+^^^,x'] dx' (3) 



Optimal coherency for high frequencies will certainly 
be achieved, if the current layer is infinitely narrow: 
j{t, x) = j{t)S{x—Xei{t)). To include more realistic cases, 
we allow in our calculations for a narrow, but finite elec- 
tron distribution: 



jit,x) ^ j{t)f{x ~ Xel{t)) 



(4) 



with variable current j{t) and position Xeiit), but 
fixed shape f{x). We take the Fourier trans- 
form of Eq. ([3]), thereby considering the re- 
tarded time, and arrive at the integral Esy{uj) — 
2ttc^^ f{ijj) J j{t) exp [— ioj + a;e((t)/c)] dt. Here, /(w) 
denotes the Fourier transform of the shape function. 
In analogy to the standard synchrotron radiation by 
a point particle, the integral can be solved with the 
method of stationary phase. Therefore we note, that 
for high u the main contributions to the integral come 
from the regions, where the phase $ = uj {t ^ Xei{t) / c) 
is approximately stationary, i.e. d^/dt w 0. However, to 
get some kind of result we need some assumption about 
the relation between the functions j{t) and Xei{t). Since 
we are dealing with the ultrarelativistic regime oq S> 1, 
it is reasonable to assume that changes in the velocity 
components are governed by changes in the direction 
of movement rather than by changes in the absolute 



velocity, which is constantly very close to the speed of 
light c. 

Let the electron motion in momentum space be given by 
iPx,Pv) = aomec{pxPy), so that j{t) = ecnePy/{aQ'^ + 
Pi +pI) ^ei = cpxl (oq ^ +v\ The derivative of 

the phase approaches zero at points where igZ ~ — c, so 
from our assumption of ultrarelativistic motion we con- 
clude = at these instants. Now, two cases have to 
be distinguished: 

1. The current changes sign at the stationary phase 
point. Then we can Taylor expand — a^t and 
Xeiit) — —vot + The integral can now be 

expressed in terms of the well-known Airy- function, 
yielding = C />) c.-2/3Ai' (^c^^/s) , 

where Ai' is the Airy function derivative and C = 
i{2Tr)^eneC^^aoai "^^^ is a complex prefactor. We 
now find the spectral envelope 
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(5) 



ai-f^, where 70 = (1 - wg) 



with LUrs ~ 2^/^ 

is the relativistic 7-factor of the electron bunch at 
the instant when the bunch moves towards the ob- 
server. As in the BGP case, the spectral envelope 
^ does not depend on all details of the electron 
bunch motion Xei , but only on its behavior close to 
the stationary points, i.e. the 7-spikes. 

2. The current does not change sign at the stationary 
phase point. Because of the assumption of highly 
relativistic motion the changes in absolute velocity 
can again be neglected compared to the changes 
in direction, and it follows that in this case the 
third derivative of Xei is zero at the stationary phase 
point. Therefore, our Taylor expansions look like 
j{t) = aot'^ and Xei{t) = —vqI + ait^ /5. This leads 



us to the spectral envelope 



/(w)cx l/Hpc^ 



-6/5 




H 2 



(6) 



with S" being the second derivative of S(a;) = 
(27r)~^ J exp [i (^xt + dt, a special case of the 

canonical swallowtail integral [U]. For the char- 
acteristic frequency ujrs we now obtain ~ 
2^/^ .^/oi^Q-^ . Because now even the derivative of 
Xei is zero at the stationary phase point, the in- 
fiuence of acceleration on the spectrum decreases 
and the characteristic frequency scaling is closer to 
the 7^-scaling for a mirror moving with constant 
velocity. 

In Fig. [4] the CSE spectra of the synchrotron radia- 
tion from the electron sheets are depicted. Comparing 
them to the 8/3-power law from the BGP-case, we no- 
tice that, because of the smaller exponents of their power 
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Figure 4: (Color online) Coherent ID synchrotron spectra for 
an infinitely thin electron layer /(tu) = f and ujrs = 100. The 
blue line corresponds to Eq. ([5| and the red line to Eq. ([6|. 
For comparison, the dashed black line denotes the BGP 8/3- 
power law. 
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Figure 5: (Color online) Normalized transverse fluid velocity 
Vy/c of the electron nanobunch. Parts of the plasma with a 
density below SOOA'^c are filtered out. 



Thus the spectral cut-ofF is determined either by cors , cor- 
responding to the relativistic 7-factor of the electrons, or 
by ujrf corresponding to the bunch width. A look at the 
motion of the electron nanobunch in the PIC simulation 
(Fig. [5| tells us that there is no change in sign of the 
transverse velocity at the stationary phase point, conse- 
quently we use Eq. ([6]). We choose w^/ = 225 Wq and 
LUrs = 800 Wo to fit the PIC spectrum, ujrf corresponds 
to a Gaussian electron bunch with a width of (5 = 10"'^ A, 
which matches reasonably well to the Sfwhm — 0.0015 A 
(see Fig. measured in the simulation, corresponding 
to a Gaussian electron bunch f{x) — exp [— (a;/^)^] with 
a width of 6 — 10~'^A and an energy of 7 ~ 10. This 
matches well with the measured electron bunch width 
(^FWHM = 0.0015 A (see figure |2]3) and the laser ampli- 
tude ao = 60, since we expect 7 to be smaller but in the 
same order of magnitude as uq. In this case utrf < ^^rs, so 
the cut-ofF is dominated by the finite bunch width. Still, 
both values are in the same order of magnitude, so that 
the factor coming from the Swallowtail-function cannot 
be neglected and actually contributes to the shape of the 
cut-ofF. The modulations that appear in Fig. |4] for fre- 
quencies around ujrs and above cannot be seen in the 
spectra, because it is suppressed by the Gauss-function 
Eq. ([7]) . The analytical synchrotron spectrum agrees ex- 
cellently with the PIC result, as the reader may verify in 
figure [T]i. 

V. PROPERTIES OF THE CSE RADIATION 

As we see in Fig. [T|;, the CSE radiation is emitted in 
the form of a single attosecond pulse whose amplitude 
is significantly higher than that of the incident pulse. 
This pulse has a FWHM duration of 0.003 laser periods, 
i.e. 9 as for a laser wavelength of 800 nm. This is very 
different from emission of the ROM harmonics, which 
need to undergo diffraction ^11) or spectral filtering j7] 
before they take on the shape of attosecond pulses. 

When we apply spectral a spectral filter in a frequency 
range (wiow, '^high) to a power-law harmonic spectrum 
with an exponent q, so that I{uj) ~ I{){ljq/uiY i the 
energy efficiency of the resulting attosecond pulse gen- 
eration process is 



law part, the CSE spectra are much flatter, around the 
100th harmonic we win more than two orders of magni- 
tude. Another intriguing property are the side maxima 
found in the spectrum ^ . This might provide an expla- 
nation for modulations that are occasionally observed in 
harmonics spectra, compare e.g. Ref. jl]. 

To compare with the PIC results, the finite size of the 
electron bunch must be taken into account. Therefore we 
assume a Gaussian density profile which leads us to 
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UJrf 



(7) 



'yatto 



^^0 

q-1 



I{uj) duj 



^0 



(8) 
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t^high 
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The scahng ^ gives rj^^^ ~ (wo/wiow)^/^ for the BGP 
spectrum with q — 8/3. For unfiltered CSE harmon- 
ics with the spectrum g = 4/3 the efficiency is close to 
^atto = 1- This means that almost the whole energy 
of the original optical cycle is concentrated in the at- 
tosecond pulse. Note that absorption is very small in 
the PIC simulations shown: it amounts to 5% in the run 
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corresponding to Fig. JTp-d and is even less in the run 
corresponding to Fig. ^-b. 

The ROM harmonics can be considered as a perturba- 
tion in the reflected signal as most of the pulse energy 
remains in the fundamental. On the contrary, the CSE 
harmonics consume most of the laser pulse energy. This 
is nicely seen in the spectral intensity of the reflected fun- 
damental for the both cases (compare Fig.[T]D and d). As 
the absorption is negligible, the energy losses at the fun- 
damental frequency can be explained solely by the energy 
transfer to high harmonics. We can roughly estimate this 
effect by /q^^^/Io^^^ « uj-^'^ dw/ uj-^/^ dw = 5. 
This value is quite close to the one from the PIC simula- 
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Further, we can estimate amplitude of the CSE at- 
tosecond pulse analytically from the spectrum. Since the 
harmonic phases are locked, for an arbitrary power law 
spectrum I{uj) oc and a spectral fllter (wiow,Whigh) 
we integrate the amplitude spectrum and obtain: 



2jl 



q-2 



^0 



f ^0 



V "^high 



(9) 



Apparently, when the harmonic spectrum is steep, i.e. 
q > 2, the radiation is dominated by the lower harmon- 
ics ojiow This is the case of the BGP spectrum q — 8/3. 
That is why one needs a spectral fllter to extract the 
attosecond pulses here. The situation changes drasti- 
cally for slowly decaying spectra with q < 2 like the CSE 
spectrum with q — 4/3. In this case, the radiation is 
dominated by the high harmonics Whigh- Even without 
any spectral flltering the radiation takes on the shape of 
an attosecond pulse. As a rule of thumb formula for the 
attosecond peak fleld of the unflltered CSE radiation we 
can write: 



77' 



CSE 



V3(r 



,1/3 



I] En 



(10) 



Using TTic — i^c/'-^o = 225, the lower of the two cut- 
off harmonic numbers used for comparison with the PIC 
spectrum in Fig. [T|l, we obtain ii'pcak = 8.8 i?o- This is 
in nice agreement with Fig. [ij;. 



VI. PARAMETRIC AL DEPENDENCE OF 
HARMONICS RADIATION 

Now, we have a look at the dependence of the har- 
monics radiation in and close to the CSE regime on the 
laser and plasma parameters. Exemplary, the laser in- 
tensity and the preplasma scale length are varied here. 
The pulse duration however will be left constantly short, 
so that we can simply focus our interest on the main op- 
tical cycle. For longer pulses, the extent of nanobunch- 
ing may vary from one optical cycle to another, which 
makes a parametrical study more difficult. We are going 
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Figure 6: Dependence of the intensity boost rj = 
max(i5^)/ ma.x{Ef) and the pulse compression F = {loot)~^ , 
where r is the FWHM width of the attosecond intensity peak 
in the reflected radiation, on ao- The laser amplitude ao is 
varied between 5 and 195 in steps of 10. Other parameters 
are the same as in Fig. |2j5. 



to examine two dimensionless key quantities: the inten- 
sity boost rj = max{E^) / max{Ef) and the pulse com- 
pression F = {ujqt)~^. It is straightforward to extract 
both magnitudes from the PIC data, and both are quite 
telling. The intensity boost is a sign of the mechanism 
of harmonics generation. If the ARP boundary condi- 
tion Eq. (jT|) is approximately valid, we must of course 
have 77 w 1. Then again, if the radiation is generated 
by nanobunches, we expect to have strongly pronounced 
attosecond peaks (compare Eq. (lOl) in the reffected ra- 
diation and therefore rj ^ 1. The pulse compression F is 
deflned as the inverse of the attosecond pulse duration. 
In the nanobunching regime, we expect it to be roughly 
proportional to rj, as the total efficiency of the attosecond 
pulse generation remains ?7atto ^ 1, compare Eq. ([s]). In 
the BGP regime, there are no attosecond pulses observed 
without spectral flltering. So the FWHM of the intensity 
peak is on the order of a quarter laser period, and we 
expect F ~ 1. 

In flgure [6] the two parameters rj and F are shown in 
dependence of ao. Except for the variation of oq, the 
parameters chosen are the same as in Figs{l]:-d, |2]d and 
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First of all we notice, that for all simulations in this 
series with ag 3> 1, we flnd r/ ^ 1. Thus, Eq. ([iJ is 
violated in all cases. Since we also notice F ^ 1 and 
F ^ T], we know, that the radiation is emitted in the shape 
of attosecond peaks with an efficiency of the order 1 . This 
indicates, that we can describe the radiation as CSE. The 
perhaps most intriguing feature of Fig. [6] is the strongly 
pronounced peak of both curves around aq = 55. We 
think that because of some very special phase matching 
between the turning point of the electron bunch and of 
the electromagnetic wave, the electron bunch experiences 
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Figure 7: Dependence of the intensity boost rj = 
max(i?^)/ ma,x{Ef) and the pulse compression F = (loot)^^ , 
on the plasma scale length in units of the laser wavelength 
L/A in the lab frame. Except for the plasma scale length, 
parameters are the same as in Fig. ^jp. The plasma ramp is 
again an exponential one oc exp{x/L). 



an unusually high compression at this parameter settings. 
This is the case that was discussed in section Hill 

Figure[7]shows the two parameters rj and T as functions 
of the plasma gradient scale length L. It is seen that both 
functions possess several local maxima. Further, j] and 
r behave similar apart from one runaway value at L = 
0.225A, where the FWHM peak duration is extremely 
short, but the intensity boost is not as high. A look at 
the actual field data tells us that in this case the foot 
of the attosecond peak is broader, consuming most of 



the energy. This deviation might e.g. be caused by a 
different, non-gaussian shape of the electron nanobunch. 

The maximum of both functions lies around L = 0.33 A, 
the parameter setting analyzed in detail before. In the 
limit of extremely small scale lengths L < O.IA, 77 and 
F become smaller, but they remain clearly bigger than 
one. Thus the reflection in this parameter range can still 
not very well be described by the ARP boundary condi- 
tion. For longer scale lengths L > 0.8A, both key values 
approach 1, so the ARP boundary condition can be ap- 
plied here. This is a possible explanation, why the BGP 
spectrum (|2| could experimentally measured at oblique 
incidence [8j. 



VII. CONCLUSIONS 

In this paper we have identified a novel mechanism 
of harmonics generation at over dense plasma surfaces, 
leading to the fiattest harmonics spectrum known so far. 
Extremely dense and narrow peaks in the electron den- 
sity, we call them nanobunches, are responsible for the 
radiation. The spectrum can be understood as a ID 
synchrotron spectrum emitted by a relativistically mov- 
ing, extremely narrow and dense electron layer. Like 
the CWE and ROM harmonics, the CSE harmonics are 
phase locked and appear in the form of attosecond pulses. 
In contrast to CWE and ROM, the attosecond pulses 
are visible immediately, therefore no energy needs to be 
wasted in spectral filtering and the energy efficiency of 
attosecond pulse generation is close to 100%. 

The work has been supported by DFG Transregio SFB 
TR18 and GraduiertenkoUeg GRK 2103. 
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